SARS-CoV-2 viral replication persists in the human lung for several weeks after onset of symptomatic severe COVID-19 and is associated with attenuated pulmonary immunity and variant-specific clinical sequalae (preprint)

The immunopathogenesis of severe COVID-19 is incompletely understood. In contradistinction to the upper respiratory tract where replicating (culturable) SARS-CoV-2 is recoverable approximately ~ 4 to 8 days after symptom onset, there is paucity of data about the frequency or duration of replicating virus in the lower respiratory tract (the human lung). We undertook lung tissue sampling (needle biopsy), within ~2 hours of death, in 42 mechanically ventilated decedents during the Beta and Delta waves. Lung biopsy cores underwent viral culture, histopathological analysis, electron microscopy, transcriptomic profiling, immunohistochemistry and cell-based flow cytometry of deconstructed tissue. 38% (16/42) of patients had culturable virus in the lung (persisting for up to 4 weeks after symptom onset). This, hitherto, undescribed bio-phenotype of lung-specific persisting viral replication was associated with an enhanced pulmonary pro-inflammatory response and variant-specific increased rates of bacterial bronchopneumonia and accelerated death. These findings question existing paradigms and suggest that in a subset of patients, concurrent, rather than sequential active viral replication continues to drive a heightened pro-inflammatory response. Our findings have potential implications for the design of therapeutic interventional strategies and clinical management of severe COVID-19 disease.

Continued Emergence and Evolution of Omicron in South Africa: New BA.4 and BA.5 lineages (preprint)

South Africa's fourth COVID-19 wave was driven predominantly by three lineages (BA.1, BA.2 and BA.3) of the SARS-CoV-2 Omicron variant of concern. We have now identified two new lineages, BA.4 and BA.5. The spike proteins of BA.4 and BA.5 are identical, and comparable to BA.2 except for the addition of 69-70del, L452R, F486V and the wild type amino acid at Q493. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure with the TaqPath COVID-19 qPCR assay. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa from the first week of April 2022 onwards. Using a multinomial logistic regression model, we estimate growth advantages for BA.4 and BA.5 of 0.08 (95% CI: 0.07 - 0.09) and 0.12 (95% CI: 0.09 - 0.15) per day respectively over BA.2 in South Africa.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance. (preprint)

Investment in Africa over the past year with regards to SARS-CoV-2 genotyping has led to a massive increase in the number of sequences, exceeding 100,000 genomes generated to track the pandemic on the continent. Our results show an increase in the number of African countries able to sequence within their own borders, coupled with a decrease in sequencing turnaround time. Findings from this genomic surveillance underscores the heterogeneous nature of the pandemic but we observe repeated dissemination of SARS-CoV-2 variants within the continent. Sustained investment for genomic surveillance in Africa is needed as the virus continues to evolve, particularly in the low vaccination landscape. These investments are very crucial for preparedness and response for future pathogen outbreaks.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function (preprint)

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Assessing the clinical severity of the Omicron variant in the Western Cape Province, South Africa, using the diagnostic PCR proxy marker of RdRp target delay to distinguish between Omicron and Delta infections – a survival analysis (preprint)

Background Emerging data suggest that SARS-CoV-2 Omicron variant of concern (VOC)is associated with reduced risk of severe disease. The extent to which this reflects a difference in the inherent virulence of Omicron, or just higher levels of population immunity, is currently not clear. Methods RdRp target delay (RTD: a difference in cycle threshold value of RdRp - E > 3.5) in the Seegene Allplex™ 2019-nCoV PCR assay is a proxy marker for the Delta VOC. The absence of this proxy marker in the period of transition to Omicron was used to identify suspected Omicron VOC infections. Cox regression was performed for the outcome of hospital admission in those who tested positive for SARS-CoV-2 on the Seegene Allplex™ assay from 1 November to 14 December 2021 in the Western Cape Province, South Africa, public sector. Vaccination status at time of diagnosis, as well as prior diagnosed infection and comorbidities, were adjusted for. Results 150 cases with RTD (proxy for Delta) and 1486 cases without RTD (proxy for Omicron) were included. Cases without RTD had a lower hazard of admission (adjusted Hazard Ratio [aHR] of 0.56, 95% confidence interval [CI] 0.34-0.91). Complete vaccination was protective of admission with an aHR of 0.45 (95%CI 0.26-0.77). Conclusion Omicron has resulted in a lower risk of hospital admission, compared to contemporaneous Delta infection in the Western Cape Province, when using the proxy marker of RTD. Under-ascertainment of reinfections with an immune escape variant like Omicron remains a challenge to accurately assessing variant virulence.

Saliva swabs are the preferred sample for Omicron detection (preprint)

The Omicron variant is characterised by more than 50 distinct mutations, the majority of which are located in the spike protein. The implications of these mutations for disease transmission, tissue tropism and diagnostic testing are still to be determined. We evaluated the relative performance of saliva and mid-turbinate swabs as RT-PCR samples for the Delta and Omicron variants. The positive percent agreement (PPA) of saliva swabs and mid-turbinate swabs to a composite standard was 71% (95% CI: 53-84%) and 100% (95% CI: 89-100%), respectively, for the Delta variant. However, for the Omicron variant saliva and mid-turbinate swabs had a 100% (95% CI: 90-100%) and 86% (95% CI: 71-94%) PPA, respectively. This finding supports ex-vivo data of altered tissue tropism from other labs for the Omicron variant. Reassessment of the diagnostic testing standard-of-care may be required as the Omicron variant becomes the dominant variant worldwide.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa (preprint)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic in southern Africa has been characterised by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, whilst the second and third waves were driven by the Beta and Delta variants respectively. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng Province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, predicted to influence antibody neutralization and spike function4. Here, we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Severity and inpatient mortality of COVID-19 pneumonia from Beta variant infection: a clinical cohort study in Cape Town, South Africa (preprint)

ABSTRACT Background The SARS-CoV-2 Beta variant, associated with immune escape and higher transmissibility, drove a more severe second COVID-19 wave in South Africa. Individual patient level characteristics and outcomes with the Beta variant are not well characterized. Methods We performed a retrospective cohort study comparing disease severity and inpatient mortality of COVID-19 pneumonia between the first and second wave periods at a referral hospital in Cape Town, South Africa. Beta variant infection was confirmed by genomic sequencing. Outcomes were analyzed with logistic regression and accelerated failure time models. Results 1,182 patients were included: 571 during the first wave period and 611 from the second wave. Beta variant accounted for 97% of infections in the second wave. There was no difference in crude in-hospital mortality between wave periods (first wave 22.2%, second wave 22.1%; p = 0.9). Time to death was decreased with higher weekly hospital admissions (16%; 95% CI, 8 to 24 for every 50-patient increase), age (18%; 95% CI, 12 to 24 for every 10-year increase) and hypertension (31%; 95% CI, 12 to 46). Corticosteroid use delayed time to death by 2-fold (95% CI, 1.5 to 3.0). Admission during the second wave decreased time to death after adjustment for other predictors, but this did not reach statistical significance (24%; 95% CI, 47 to -2). There was no effect of HIV on survival. Conclusions There was a trend towards earlier mortality during the second COVID-19 wave driven by the Beta variant, suggesting a possible biological basis. Use of oral prednisone was strongly protective. Key points In Cape Town, South Africa, the second wave of COVID-19, dominated by the Beta variant, was associated with decreased time to inpatient death after adjustment for age, comorbidities, steroid use, and admission numbers. Use of oral prednisone was strongly protective.

Higher mortality associated with the SARS-CoV-2 Delta variant in the Western Cape, South Africa, using RdRp target delay as a proxy (preprint)

A novel proxy for the Delta variant, RNA-dependent RNA polymerase target delay in the Seegene Allplex™ 2019-nCoV PCR assay, was associated with higher mortality (adjusted Odds Ratio 1.45 [95%CI 1.13-1.86]), compared to presumptive Beta infection, in the Western Cape, South Africa (April-July 2021). Prior diagnosed infection and vaccination were protective.

Reduced amplification efficiency of the RNA-dependent-RNA-polymerase (RdRp) target enables tracking of the Delta SARS-CoV-2 variant using routine diagnostic tests (preprint)

Routine SARS-CoV-2 surveillance in the Western Cape region of South Africa (January-August 2021) found a reduced PCR amplification efficiency of the RdRp gene target of the Seegene, Allplex 2019-nCoV diagnostic assay when detecting the Delta variant. We propose that this can be used as a surrogate for variant detection.

Rapid replacement of the Beta variant by the Delta variant in South Africa (preprint)

The Beta variant of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in South Africa in late 2020 and rapidly became the dominant variant, causing over 95% of infections in the country during and after the second epidemic wave. Here we show rapid replacement of the Beta variant by the Delta variant, a highly transmissible variant of concern (VOC) that emerged in India and subsequently spread around the world. The Delta variant was imported to South Africa primarily from India, spread rapidly in large monophyletic clusters to all provinces, and became dominant within three months of introduction. This was associated with a resurgence in community transmission, leading to a third wave which was associated with a high number of deaths. We estimated a growth advantage for the Delta variant in South Africa of 0.089 (95% confidence interval [CI] 0.084-0.093) per day which corresponds to a transmission advantage of 46% (95% CI 44-48) compared to the Beta variant. These data provide additional support for the increased transmissibility of the Delta variant relative to other VOC and highlight how dynamic shifts in the distribution of variants contribute to the ongoing public health threat.

The continuous evolution of SARS-CoV-2 in South Africa: a new lineage with rapid accumulation of mutations of concern and global detection (preprint)

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The emergence of C.1.2, associated with a high substitution rate, includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 VOC/VOIs. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta showed high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

Prior infection with SARS-CoV-2 boosts and broadens Ad26.COV2.S immunogenicity in a variant dependent manner (preprint)

The Johnson and Johnson Ad26.COV2.S single dose vaccine, designed as an emergency response to the pandemic, represents an attractive option for the scale-up of COVID-19 vaccination in resource-limited countries. We examined the effect of prior infection with ancestral (D614G) or Beta variants on Ad26.COV2.S immunogenicity approximately 28 days post-vaccination. We compared healthcare workers who were SARS-CoV-2 naive (n=20), to those infected during the first wave prior to the emergence of Beta (n=20), and those infected in the second wave (n=20), when Beta was the dominant variant. We demonstrate that a priming exposure from infection significantly increased the magnitude of spike binding antibodies, neutralizing antibodies and antibody-dependent cellular cytotoxicity activity (ADCC) against D614G, Beta and Delta variants. The magnitude of antibody boosting was similar in both waves, despite the longer time interval between wave 1 infection and vaccination (7 months), compared to wave 2 (2 months). ADCC and binding cross-reactivity was similar in both waves. However, neutralization cross-reactivity varied by wave, showing that the antibody repertoire was shaped by the spike sequence of the infecting variant. Robust CD4 and CD8 T cell responses to spike of similar or higher magnitude as those elicited by infection were induced after vaccination. In contrast to antibody responses, prior infection was not required for the generation of high magnitude T cell responses, and T cell recognition of the Beta variant was fully preserved. Therefore, Ad26.COV2.S vaccination following prior infection, even >6 months previously, may result in substantially enhanced protection against COVID-19, of particular relevance in settings of high SARS-CoV-2 seroprevalence. Furthermore, the dominant impact of the infecting variant on neutralization breadth after vaccination has important implications for the design of second-generation vaccines based on variants of concern.

Loss of recognition of SARS-CoV-2 B.1.351 variant spike epitopes but overall preservation of T cell immunity (preprint)

SARS-CoV-2 variants have emerged that escape neutralization and potentially impact vaccine efficacy. T cell responses play a role in protection from reinfection and severe disease, but the potential for spike mutations to affect T cell immunity is poorly studied. We assessed both neutralizing antibody and T cell responses in 44 South African COVID-19 patients infected either with B.1.351, now dominant in South Africa, or infected prior to its emergence (first wave), to provide an overall measure of immune evasion. We show for the first time that robust spike-specific CD4 and CD8 T cell responses were detectable in B.1.351-infected patients, similar to first wave patients. Using peptides spanning only the B.1.351 mutated regions, we identified CD4 T cell responses targeting the wild type peptides in 12/22 (54.5%) first wave patients, all of whom failed to recognize corresponding B.1.351-mutated peptides (p=0.0005). However, responses to the mutated regions formed only a small proportion (15.7%) of the overall CD4 response, and few patients (3/44) mounted CD8 responses that targeted the mutated regions. First wave patients showed a 12.7 fold reduction in plasma neutralization of B.1.351. This study shows that despite loss of recognition of immunodominant CD4 epitope(s), overall CD4 and CD8 T cell responses to B.1.351 are preserved. These observations may explain why, despite substantial loss of neutralizing antibody activity against B.1.351, several vaccines have retained the ability to protect against severe COVID-19 disease.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa (preprint)

The progression of the SARS-CoV-2 pandemic in Africa has so far been heterogeneous and the full impact is not yet well understood. Here, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations, predominantly from Europe, which diminished following the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1 and C.1.1. Although distorted by low sampling numbers and blind-spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a breeding ground for new variants.

Improved performance of saliva for the detection of the B.1.351 variant in South Africa (preprint)

Assessment of the unknown performance of saliva for the detection of the novel SARS-CoV-2 variant of concern (VOC) B.1.351 (501Y.V2) lineage is essential as saliva has been shown to be an equivalent, less invasive and a less costly alternative to nasopharyngeal swabs for the molecular detection of SARS-CoV-2 infection in pre-variant studies. Between 1st August 2020 and 16th January 2021, we enrolled 410 eligible ambulatory participants who presented to Groote Schuur Hospital (GSH) in Cape Town, South Africa for SARS-CoV-2 testing. Of these, 300 were enrolled prior to, and 110 after, the initial detection and replacement of wild-type by the B.1.351 variant. All participants provided a supervised self-collected mid-turbinate (MT) and saliva (SA) swab, in addition to the standard HCW collected NP swab which were all tested by RT-PCR in an accredited diagnostic laboratory. Positive percent agreement to NP swab for SA swabs pre- and post-variant were 51.5% and 72.5% respectively while these values for MT swabs were 75.8% and 77.5%. The negative percent agreement for all swab types during all periods was >98%. The basis for this marked improvement of SA swabs as a diagnostic sample for B.1.351 virus is still being investigated.

SARS-CoV-2 501Y.V2 (B.1.351) elicits cross-reactive neutralizing antibodies (preprint)

Neutralization escape by SARS-CoV-2 variants, as has been observed in the 501Y.V2 (B.1.351) variant, has impacted the efficacy of first generation COVID-19 vaccines. Here, the antibody response to the 501Y.V2 variant was examined in a cohort of patients hospitalized with COVID-19 in early 2021 - when over 90% of infections in South Africa were attributed to 501Y.V2. Robust binding and neutralizing antibody titers to the 501Y.V2 variant were detected and these binding antibodies showed high levels of cross-reactivity for the original variant, from the first wave. In contrast to an earlier study where sera from individuals infected with the original variant showed dramatically reduced potency against 501Y.V2, sera from 501Y.V2-infected patients maintained good cross-reactivity against viruses from the first wave. Furthermore, sera from 501Y.V2-infected patients also neutralized the 501Y.V3 (P.1) variant first described in Brazil, and now circulating globally. Collectively these data suggest that the antibody response in patients infected with 501Y.V2 has a broad specificity and that vaccines designed with the 501Y.V2 sequence may elicit more cross-reactive responses.

Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa (preprint)

Continued uncontrolled transmission of the severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in many parts of the world is creating the conditions for significant virus evolution. Here, we describe a new SARS-CoV-2 lineage (501Y.V2) characterised by eight lineage-defining mutations in the spike protein, including three at important residues in the receptor-binding domain (K417N, E484K and N501Y) that may have functional significance. This lineage emerged in South Africa after the first epidemic wave in a severely affected metropolitan area, Nelson Mandela Bay, located on the coast of the Eastern Cape Province. This lineage spread rapidly, becoming within weeks the dominant lineage in the Eastern Cape and Western Cape Provinces. Whilst the full significance of the mutations is yet to be determined, the genomic data, showing the rapid displacement of other lineages, suggest that this lineage may be associated with increased transmissibility.

Major new lineages of SARS-CoV-2 emerge and spread in South Africa during lockdown. (preprint)

In March 2020, the first cases of COVID-19 were reported in South Africa. The epidemic spread very fast despite an early and extreme lockdown and infected over 600,000 people, by far the highest number of infections in an African country. To rapidly understand the spread of SARS-CoV-2 in South Africa, we formed the Network for Genomics Surveillance in South Africa (NGS-SA). Here, we analyze 1,365 high quality whole genomes and identify 16 new lineages of SARS-CoV-2. Most of these unique lineages have mutations that are found hardly anywhere else in the world. We also show that three lineages spread widely in South Africa and contributed to ~42% of all of the infections in the country. This included the first identified C lineage of SARS-CoV-2, C.1, which has 16 mutations as compared with the original Wuhan sequence. C.1 was the most geographically widespread lineage in South Africa, causing infections in multiple provinces and in all of the eleven districts in KwaZulu-Natal (KZN), the most sampled province. Interestingly, the first South-African specific lineage, B.1.106, which was identified in April 2020, became extinct after nosocomial outbreaks were controlled. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify and control the spread of SARS-CoV-2.